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A tribological system consisting of two contacting blocks has been considered. One of them 
is arranged between two springs, the other is driven periodically. The kinetics of the system 
has been studied in the boundary friction mode, when an ultrathin lubricant film is contained 
between the atomically smooth surfaces. In order to describe the film state, the expression for 
the free energy density is used in the form of an expansion in a power series in the order 
parameter, the latter being reduced to the shear modulus of a lubricant. The stick-slip mode 
is shown to be realized in a wide range of parameters, being a result of the periodic first-order 
phase transitions between kinetic friction regimes. The behavior of the system governed by 
internal and external parameters has been predicted. 
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1. Introduction 



Owing to the rapid development of high-precision 
experimental techniques aimed at researching the 
nanosystems, the processes of friction in the case 
where the thickness of a lubricant between rubbing 
surfaces equals several atomic diameters have been 
intensively studied recently (H-Q- I n tribology, this 
friction mode was called "boundary friction". It is 
often realized in ordinary mechanisms, because the 
rubbing surfaces contact with each other by means 
of surface irregularities or inhomogeneities [l|,|4( . The 
boundary friction differs essentially from the hydro- 
dynamic mode, when the friction force is proportional 
to a power function of the velocity. Note that the 
ultrathin film of a lubricant does not form conven- 
tional, thermodynamically equilibrium phases, solid 
and liquid ones. Instead, we have liquid- and solid- 
like states, which are kinetic friction modes, and there 
can be several of them [H,[6j]. This occurs because the 
symmetry of a lubricant state is substantially affected 
by friction surfaces, this fact being not of importance 
for bulk lubricants. In the course of friction, the phase 
transitions of both the first and second orders can 
take place between the stationary states ■ These 
transitions often comprise the origin of the stick-slip 
motion mode for contacting surfaces [5|,|6|,|9j. 
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In order to describe the boundary mode of fric- 
tion and nano-contact phenomena, phenomenological 
models are widely used (7l4loj|. In particular, a model 
was developed in which the lubricant melting is 
driven by the thermodynamic and shear mechanisms. 
In the framework of this model, the influence of ad- 
ditive fluctuations of principal quantities was stud- 
ied [ijj], and their presence in the system was demon- 
strated to result in the emergence of new stationary 
states and new kinetic friction modes [l3, 14 1, which 
are not essential for bulk systems. The origin of the 
hysteretic behavior, which was observed experimen- 
tally [l,[lM[iB|) was elucidated in work [l7|. The in- 
dicated model also made it possible to describe the 
periodic stick-slip mode of motion 3, lj| . 

In works 0, H3; a thermodynamic scenario of 
boundary friction was proposed, which is based on 
the phase transition theory developed by Landau [2l[ . 
This model takes into account that the ultrathin 
film of a lubricant can melt and stay in a liquid- 
like disordered state both owing to the ordinary ther- 
modynamic melting and as a result of overcoming 
the fluidity threshold by the shear stress component 
("shear melting"). The influence of those factors was 
also studied in work |8j, in which the excess vol- 
ume arising owing to the lubricant stochas- 
tization at its melting was selected as the order pa- 
rameter. As the excess volume increases, the shear 



ISSN 2071-0186. Ukr. J. Phys. 2013. Vol. 58, No. 1 



91 



LA. Lyashenko, A.M. Zaskoka 



K 









M 



K 



Fig. 1. Diagram of a tribological system 

modulus decreases [H, which results in the melting. 
In works @,[2(| , the shear modulus itself was selected 
as the order parameter, which acquires zero values 
in the liquid- like phase. However, in works [7LI20I. 
the melting is described as a continuous phase tran- 
sition of the second order, whereas jump-like phase 
transitions of the first order are often observed in the 
boundary friction mode , which are responsible 

for the stick-slip motion [5|,[6j]. 

This work is aimed at describing the first-order 
phase transition in the framework of the model de- 
veloped in works [7,20] and at studying the behavior 
of tribological systems on the basis of the indicated 
modification. The model proposed does not make al- 
lowance for specific types of lubricants, because its 
task consists in describing the origins of phenomena 
that take place at the boundary friction. For specific 
types of lubricants and friction surfaces, the model 
should be modified. To some extent, it can be made 
by choosing the numerical values of coefficients in 
the series expansions of the free energy, relaxation 
times, and so on. The model describes only homoge- 
neous lubricants composed of non-polar quasispher- 
ical molecules [5, 6]. One of the reasons is the fact 
that we study a situation where the elastic stresses 
acquire zero values in the liquid-like state, i.e. the 
melting gives rise to the total disordering of lubricant 
molecules, which does not takes place in thin lubri- 
cant films consisting of polymer molecules. Another 
reason is the fact that the obtained time dependences 
of the friction force and stresses are strictly periodic, 
which is also observed for quasispherical molecules 
only [113. 

2. Tribological System 

The fabrication of atomically smooth surfaces with 
large dimensions is associated with considerable tech- 
nological difficulties. Therefore, to measure the dy- 
namic parameters of ultrathin lubricant films between 
such surfaces, the surfaces characterized by small di- 
mensions and pasted on spherical or cylindrical sur- 
faces that rub each other are used. This scheme was 
applied while designing the surface force apparatus 



(SFA) dHUH. Two types of SFA-Mk II and 
Mk III— were described in review In the latter, 
the system to control distances between rubbing sur- 
faces was improved. The device allows the shape of 
surfaces to be determined, as well as the distance be- 
tween them to within an accuracy of 1 A. The contact 
area between surfaces is measured with an accuracy of 
±5%, the normal and shear components of operating 
forces to within ±1%, and the magnitude of applied 
loading to within ±5%. 

One of the rubbing surfaces in the SFA is fixed, 
and the other is driven to move periodically. In the 
course of motion, the shear stresses and the effec- 
tive viscosity of a lubricant are measured, the lubri- 
cant structure is determined, and so forth. In this 
work, we consider a simplified mechanical analog of 
the SFA exhibited in Fig. 1. Two springs character- 
ized by the stiffness constant K are connected with 
a block of mass M mounted on rollers. The rolling 
friction for the latter is neglected below. On the in- 
dicated block, another block is arranged, which is 
brought into a periodic motion by applying an exter- 
nal force. Provided that the surfaces of two blocks 
interact with each other, the motion of the upper 
block stimulates the motion of the lower one. The 
trajectory of the lower block substantially depends 
on the friction mode established in the system. A 
similar tribological system was experimentally stud- 
ied in works 2g,22|. Note that, in contrast to the 
SFA design, now both blocks are mobile, which en- 
ables the time dependences of block coordinates and 
velocities to be registered, and, by analyzing them, 
the rheological and tribological characteristics of the 
system to be determined. 

Let X and V — X be the coordinate and the veloc- 
ity, respectively, of the upper block, whereas x and 
v = x denote the corresponding quantities for the 
lower block. Let us consider the case where the up- 
per block moves according to the cyclic law, 



X = X m coscji, 



V 



-X m u sin ut, 



(1) 
(2) 



where X m is the amplitude, and oj is the cyclic fre- 
quency. We write down the equation of motion for 
the lower block in the form [26] 



Mx + 2Kx - F = 0, 



(3) 



where F is the friction force that arises between the 
blocks at their relative motion. From the last expres- 
sion, it follows that the character of motion in the 
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system essentially depends on the friction mode and 
the lubricant properties, because they determine the 
force F. 

The friction force is determined in a standard way, 



F = a A, 



(4) 



where a is the shear stress that arises in the lubricant, 
and A is the contact area between rubbing surfaces. 

In the boundary friction mode, the elastic, a e \, and 
viscous (dissipative), a v , stresses arise in the lubri- 
cant layer [7l.l8l.ll6j|. As a rule, the melting is accom- 
panied by a reduction of the elastic stress component, 
whereas the viscous one grows owing to an increase 
of the relative shear velocity between rubbing sur- 
faces [l6j . Hence, the total stress is determined by 
the sum of indicated components, 



Cel 



(5) 



The viscous stresses in the lubricant layer are deter- 
mined by the empirical formula [28, 29] 



VeffiV - v) 



(6) 



where the effective viscosity of the lubricant, r] e g (it 
depends on plenty of factors and is determined exper- 
imentally) is introduced into consideration, as well as 
the relative velocity of surface motion, V — v. 

As a result, in the case of boundary friction, poly- 
mer solutions or melts are applied as lubricants. The 
necessity of such an application is caused by the fact 
that the friction surfaces are small in dimensions, and 
the lubricant film between them must not be squeezed 
out under the influence of large tribological loadings. 
Such lubricants are non-Newtonian fluids, the viscos- 
ity of which depends not only on the temperature, but 
also on the velocity gradient. However, the applica- 
tion of the SFA allows the behavior of a wide class 
of lubricants to be examined in the boundary friction 
mode, because the rubbing surfaces in these experi- 
ments are completely imbedded into a vessel with a 
liquid to study, so that the latter is not squeezed out 
from the gap between the surfaces during their mo- 
tion However, we should note that even ordinary 
water, when being used as a boundary lubricant, can 
behave as a non-Newtonian fluid, because, owing to 
its interaction with the surfaces, it can create spa- 
tially ordered structures in the course of motion. 

The non-Newtonian fluids are divided into two 
classes: pseudoplastic fluids, the viscosity of which 



decreases with the growth of the strain rate (e.g., 
these are polymer solutions and melts) and dilatant 
ones, the viscosity of which increases as e grows (e.g., 
suspensions of solid particles) . For both situations to 
be taken into account, let us use a simple power-law 
approximation [2Sl . [2Q| | 



(7) 



Here, we introduced the proportionality coefficient k 
(its dimension is Pa • s 7+1 ) and the dimensionless in- 
dex 7 (for pseudoplastic fluids, 7 < 0; dilatant ones 
are characterized by the index 7 > 0; and 7 = in 
the case of Newtonian fluids). 

The strain rate is determined through the relative 
velocity of motion and the lubricant thickness h [28| , 



V-v 



(8) 



Taking Eqs. ([7]) and ([5]) into account, the expression 
for viscous stresses (Eq. ([5])) looks like 



= k 



V-v 



7+1 



(9) 



Note that, according to Eq. ([9]), viscous stresses are 
available in both the liquid- and solid-like states of 
a lubricant. The presence of viscous (dissipative) 
stresses in bo th p hases was indicated in the exper- 
imental work [16]. However, if the lubricant is in a 
solid-like state, viscous stresses are low, because, in 
accordance with Eq. ^ , they are proportional to the 
relative shear velocity, V — v, which is low in this case. 

Substituting Eqs. (JSj and © into Eq. g]), we ob- 
tain the final expression for the friction force 30l-33|. 



F 



(Tci + k sgn(V - v) 



V 



7+1 



-4, 



where the function 

'1, V>v, 



sgn(V — v) = 



-1, V < v 



(10) 



(11) 



takes into account the direction of force action. The 
first term in Eq. (|10|) describes the elastic component 
of the friction force, and the second is responsible 
for the viscous one, which grows with the velocity. 
Hence, the friction force depends on the velocity of 
the lower block, v, and the elastic stresses, cr e i, that 
arise in a lubricant. 
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3. Thermodynamic Model 

In the homogeneous case, the free energy density for 
an ultrathin lubricant layer looks like [7 , 2^, 3(J 31 1 



4- (m m \ 2 i a 2 2 ^ 3 , C 4 

f = a{T - T c )tp + ~tp £ el - -tp* + -tp , 



(12) 



where T is the lubricant temperature; T c is the crit- 
ical temperature; e e i is the shear component of the 
elastic strain; a, a, 6, and c are positive constants, 
and ip is the order parameter (the amplitude of the 
periodic component in the microscopic medium den- 
sity function @,0|). The parameter ip equals zero in 
the liquid-like phase and acquires nonzero values in 
the solid- like one. In comparison with works [7l.[20|. 
potential (fT2"|) takes additionally the third-order term 
into account. This form of expansion is used to de- 
scribe phase transitions of the first order [21 , 34 1 . In 
the second term in Eq. (fT2"l) . we also introduced the 
factor a, which allows us to vary the contribution of 
the elastic energy to the potential. 

The elastic stresses that arise in the lubricant layer, 
according to Eq. (fT2|) . are determined as the deriva- 
tive (T e i = df /de c i, so that 



Oci = atp s e \. 



(13) 



Therefore, after the coefficient a has been introduced 
into expansion (|12p . the shear modulus is determined 
as follows: 



/! = dip . 



(14) 




Fig. 2. Dependences of the free energy density / (see Eq. 1121 1) 
on the dimensionless order parameter tp for various temper- 
atures T = 265, 286, and 310 K (curves 1 to 3, respec- 
tively). The calculation parameters are a = 0.95 J/(K ■ m 3 ), 
T c = 290 K, a = 4 X 10 12 Pa, b = 230 J/m 3 , c = 850 J/m 3 , 
and the shear strain e c i = 2.1 X 10 -6 



Depending on the value of parameter a, it can acquire 
both small and large values at \ip\ < 1. Note that, in 
the boundary friction mode, the shear modulus can 
be several orders of magnitude larger than that in 
the hydrodynamic mode for the same lubricant. As 
a result, if the critical temperature T or critical elas- 
tic shear stress <7 i become exceeded in the course 
of friction, the lubricant does not melt completely; 
instead, a domain structure with regions of liquid- 
assisted and dry friction is created. For this situation 
to be studied, Eq. (fT2j) must include gradient terms, 
which considerably complicates the subsequent con- 
sideration. However, the examination of such spatial 
structures comprises a separate problem, which is not 
the purpose of this work. Therefore, the gradient 
terms are excluded from Eq. (1121) . which corresponds 
to the consideration of a lubricant behavior in the 
framework of the one-domain model with a homoge- 
neous structure. 

According to the principle of minimum energy, the 
system tends to occupy a stationary state, which cor- 
responds to the minimum of the free energy f(tp) (see 
Eq. lfT2|) ). irrespective of its initial conditions. Since 
the parameter ip is the amplitude of the periodic com- 
ponent in the microscopic medium density function, 
we consider below only the physical range of values 
tp > 0. Let us introduce the function 



B{s eh T) = asly 



2a(T-T c ). 



(15) 



The analysis of expression (|12|) for the free energy al- 
lows the following situations to be distinguished. Pro- 
vided that the condition B(e e \, T) < is obeyed, the 
maximum of potential (|12p at <p = and its min- 
imum at tp > are realized (curve 1 in Fig. 2). 
In this case, the lubricant is solid-like, because the 
shear modulus /i > 0. In the intermediate interval 
< B(e c \,T) < 6 2 /(4c), the maximum of the poten- 
tial at tp = transforms into a minimum and, ad- 
ditionally, there emerges a maximum that separates 
the zero and nonzero minima (curve 2 in Fig. 2; it is 
also shown scaled- up in the inset). In this case, the 
state of a lubricant depends on the initial conditions, 
and the lubricant can be in either solid- or liquid-like 
state. In the latter case, B(e c \,T) > b 2 /(Ac), and a 
single minimum of the potential at tp = is realized 
(curve 3 in Fig. 2), which, according to Eq. (fl4|) . cor- 
responds to the zero value of lubricant shear modulus 
and its liquid-like structure. 

The stationary values of order parameter tp are 
determined as the roots of the equation df /dtp — 
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namely, 




- -et 



2a{T - T c ) 



(16) 



The root is related to the unstable stationary 
state, because it corresponds to the maximum of po- 
tential (TT2"|) . The stable state, which corresponds to 
the potential minimum, is given by the root 9?+ . Be- 
sides roots (TT6"]) , the stationary solution ipo — always 
exists, which corresponds to the extremum of poten- 
tial p2p with the order parameter equal to zero; it can 
be either a maximum or minimum of the potential. 
According to Eq. (|16[) . the lubricant melts if either 
the temperature T elevates or the shear component 
of the elastic strain, e e \, grows. Thus, the model con- 
cerned makes allowance for both thermodynamic and 
shear meltings. 

As was already indicated above, at small values 
of temperature T and strain e e \, when the function 
B{e \,T) < 0, the lubricant is solid-like, because, in 
accordance with Eq. (TT6"]) , a stationary value of pa- 
rameter tp different from zero is realized, and, accord- 
ing to Eq. HH), the shear modulus [i is also non-zero. 
In this case, the potential has a single minimum at 
tp > 0. If the temperature T exceeds the critical value 



Too — T c — 



2a 



8ac 



(17) 



the order parameter vanishes in a jump-like manner, 
when the lubricant passes into the liquid-like state, in 
which the potential /(</?) has a single minimum at ip = 



[30,31]. If, after this transition, the temperature T 
falls down further, the lubricant solidifies following 
the mechanism of first-order phase transformation at 
a lower temperature, 



T> = T r 



2a 



'el 1 



(18) 



and the parameter tp becomes non-zero again. In the 
intermediate temperature region, T~ < T < T c0l the 
potential is characterized by two minima at positive 
tp. Hence, the dependence <p(T) has a hysteresis char- 
acter [30II3H and corresponds to the phase transition 
of the first order. Expression (TTB")) elucidates the phys- 
ical meaning of the critical temperature T c ; namely, 
it is the temperature of lubricant solidification at zero 
strains, when only the mechanism of thermodynamic 
melting is active in the system. 

From expression (flTf . it follows that the lubricant 
melts not only at the temperature elevation, but also 



if it is subjected to an external mechanical action, 
when the elastic strain component exceeds the critical 
value 



£el,cO 



2a(T c - T) 



4ac' 



(19) 



Using formula (|18[) . we can determine the elastic de- 
formation £ c i, at which the lubricant solidifies, 



.0 



2a(T c - T) 



(20) 



Note that, according to relation (flU)) . the melting can 
occur even at the zero temperature, T = 0, if the 
strain exceeds the critical value. At the zero strain, 
i.e. at e c i = 0, the lubricant melts, when its temper- 
ature exceeds the critical value T c q (see Eq. (|17[l). 

As a rule, it is the relative shear velocity between 
friction surfaces rather than the shear strain compo- 
nent e c i that is registered in experiments [5,6j. There- 
fore, for our research to go further, it is necessary to 
obtain a relation between those two quantities. Let us 
take advantage of the Debye approximation, accord- 
ing to which the elastic strain component e e \ arises in 
the lubricant layer, when the latter flows plastically 
at the velocity [Tj 



Eel 
7V 



(21) 



where r 6 is the Maxwell relaxation time for internal 
stresses. The total strain in the layer is determined 
as the sum of elastic, e c i, and plastic, e p i, compo- 
nents 0,111] 



S — £cl + £pl- 



(22) 



Combining relations (JSJ, ([2"Tj) . and (|2"2"|) , we obtain 
the kinetic equation for the evolution of the elastic 
component of the shear strain 0, 3^, 31. 33 1 : 



-Eel 



(V - v)n 



(23) 



Boundary friction experiments testify that the re- 
laxation time of the elastic strain is very short, as a 
rule. This quantity can be estimated from the relation 
t £ sa a/c ~ 10~ 12 s, where a ~ 1 nm is the lattice con- 
stant or the intermolecular distance, and c ~ 10 3 m/s 
is the sound velocity [111 ]. However, in the boundary 
mode, the strain relaxation time, r e , can differ by sev- 
eral orders of magnitude Bearing in mind that 
the value of strain relaxation time t e is small, below 
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we use the adiabatic approximation r £ e c i « [35], 
which allows us to determine the strain by its sta- 
tionary value 



(V 



(24) 



rather than by Eq. 

In the general case, the free energy (fl2"j) depends on 
the lubricant layer thickness h [20]. Note that, in the 
framework of our model, the second term in expres- 
sion (JTSJ) is proportional to the square of the elastic 



strain, e 2 v In accordance with relation (|24|) . the sta- 
tionary elastic strain increases with a reduction of the 
lubricant thickness h. Therefore, in the limiting case 
of a very thin layer (h — > 0), the strain e c i — > oo. 
In this case, the second term in expansion (|12l) dom- 
inates, and the stationary value of order parameter- 
equals zero, so that the lubricant is liquid-like, as in 
work [2(| . A detailed study of the influence of the lu- 
bricant layer thickness on friction modes was carried 
out in works 36, 371 . 



4. Kinetics of Melting 

The changes in the lubricant temperature T and the 
strain e e i induce variations of the order parameter 
f, which governs the free energy / (see Eq. (j!2|) ) 
in accordance with the power-law expansion of the 
latter [2lfl. The stabilization time for a new sta- 
tionary value f + (see Eq. ([To]) ) is determined by 
the generalized thermodynamic force —df /dtp. If 
if w ip + . this force is small, and the relaxation pro- 
cess is described by the Landau-Khalatnikov linear 
kinetic equation [38| 



Of 



(25) 



where the kinetic coefficient <5 characterizes the iner- 
tial properties of the system. After substituting en- 
ergy (TT2"]) into Eq. ([23)) , we obtain the equation in the 
explicit form, 

f = -5 (2a(T - T c )f + afs 2 cl - bf 2 + cf 3 ) + 

(26) 

Equation (|2"6"]l includes a term responsible for ad- 
ditive fluctuations with a low intensity 33,31|. Their 
intensity is selected to be so low that they do not 
affect the deterministic behavior of the system. How- 
ever, their introduction is necessary, because, at sub- 
sequent numerical calculations, the root tp = of 



Eq. (j26p corresponding to the maximum of the poten- 
tial f(f), i.e. to the unstable stationary state, turns 
out stable according to the structure of the equation. 
In this situation, the introduction of £(t) stimulates 
the system to transit from the unstable state into 
a stable one, which corresponds to the energy mini- 
mum. Hence, fluctuations are taken into considera- 
tion by means of the features in subsequent numerical 
calculations. 

The dynamic characteristics of any tribological sys- 
tem are governed by its properties in whole. For 
instance, in the geometry illustrated in Fig. 1, the 
behavior of the system substantially depends on the 
stiffness constant, K, of the spring and the mass of 
the lower block, M. In contrast to the case of motion 
with constant elastic strains, this tribological system 
can reveal the stick-slip mode of motion in the course 
of friction [5, 6, 9, 18]. The indicated mode is estab- 
lished because the lubricant periodically melts and 
solidifies in the course of motion, which leads to the 
oscillatory character of friction force F. To calculate 
the evolution of this system in time, we need to solve 
the system of kinetic equations ([3]) and (|26|) numeri- 
cally, determining the friction force F from Eqs. (flQ|) 
and ([TT|) . the elastic stresses a c \ from Eq. (fiU)) . and 
the strain e e \ from relation (|24D . In so doing, we have 
to take into account the relation x — v, as well as 
definitions (1) and (2). 

While solving the differential equations numeri- 
cally, we used the Euler-Cramer method with the 
time increment At = 10~ 10 s. The initial conditions 
fa = xq = vq = were chosen. The result obtained 
is shown in Fig. 3. The dashed curve in the upper 
panel corresponds to the time dependence of the up- 
per block coordinate X(t) (see Eq. (1)), and the solid 
one to that of the lower block, x(t), which is more 
complicated. The figure also exhibits the time depen- 
dences for the block velocities, elastic shear stresses 
<j e \ (see Eq. (]13[) ) that arise in the lubricant, and total 
friction force F (see Eq. (flQ|)). Let us examine these 
dependences in more details. 

At the initial time moment, t = 0, the blocks are 
motionless, and the lubricant is in the solid-like state, 
because the dependences are plotted for the lubricant 
temperature T lower than the critical one, Tj? (see 
Eq. (|18|)). and e c \ — at rest. At t = 0, the upper 
block starts to move, and, at t > 0, its velocity grows 
in accordance with Eq. (2). Since the lubricant is in 
the solid-like state, the friction force F possesses both 
the viscous and elastic components, and the lower 
block moves together with the upper one. However, 
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Fig. 3. Dependences of the coordinates X and x, the veloc- 
ities V and ti, the elastic stresses a e \ (see Eq. ([13}), and the 
friction force F (see Eq. 1101 1 on the time t for the same pa- 
rameters as in Fig. 2 and h = 10 — 9 m, t s = 10~ 8 s, 7 = —2/3, 
A = 0.6 x 10" 9 m 2 , fe = 5 x 10 4 Pa ■ s 1 / 3 , 5 = 100 m 3 /(J ■ s), 
T = 200 K, X m = 0.8 x 10" 6 m, u = 10 rad/s, Af = 0.4 kg, 
and is* = 3000 N/m. The dashed curves correspond to the co- 
ordinate X(t) and the velocity V(t) of the upper block, and 
the solid ones to the coordinate x(t) and the velocity v(t) of 
the lower block 



in the course of motion, the absolute value of elastic 
force 2Kx, which impedes the lower block to move, 
grows and, as a result, the velocity v does not increase 
so sharply as the velocity V does. Hence, the relative 
shear velocity between block surfaces, V—v, increases 
in time and, in accordance with Eq. (|24[) . the elastic 
strain e e i also grows. At a certain time moment, the 
condition e c i > e c i,co ( see Eq- (H2)) becomes satisfied, 
and the lubricant begins to melt following the "shear 
melting" mechanism. The friction force substantially 
decreases at that, because the stresses vanish, and 
the lower block can slide for a considerable distance, 
being driven by the elastic force from the compressed 
and stretched springs. Therefore, the relative shear 




Fig. 4- Dependences of the relative displacement, X — x, 
and the velocity, V — v, of blocks on the time t for the same 
parameters as in Fig. 3 

velocity diminishes, and, when the condition e e i < 
e ei c ( see Eq. is obeyed, the lubricant solidifies 

again. The considered process repeats periodically. 

In addition, Fig. 4 demonstrates the time depen- 
dences of the relative block displacement and the 
velocity. At the time moments, when the surfaces 
"stick" to each other, their relative displacement X—x 
remains constant, and the relative shear velocity V—v 
is close to zero (in this case, the dependences V(t) 
and v(t) in Fig. 3 visually coincide). Hence, the pe- 
riodic stick-slip mode of motion takes place, which 
is also ty pica l of dry friction, when no lubricant is 
used P,L2, 39]. For the chosen parameter values, the 
blocks "stick" to each other four times during a com- 
plete period of parameter changes: two times in each 
direction of motion, with the obtained dependences 
being symmetric with respect to the motion direc- 
tion. However, a number of different situations can 
be realized depending on the system parameters. 

The phase portraits of the system calculated at the 
same parameters as in Fig. 3 and various values of 
cyclic frequency u) are depicted in Fig. 5. The ki- 
netic dependences in Fig. 3 completely correspond to 
the phase portrait in Fig. 5, a, because they were cal- 
culated for the same frequency u> value. It is im- 
portant to emphasize the fact that the phase por- 
traits in Figs. 5, a and c are symmetric with respect 
to the coordinate origin, whereas the phase portraits 
in Fig. 5,6 and d illustrate the situation where the 
motion of the upper friction surface in one direction 
affects differently the motion of the lower block in 
comparison with its motion in the opposite direction. 
Hence, the system reveals memory effects, which were 
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Fig. 5. Phase portraits of the system at the same parameters as in Fig. 3 and various values of the cyclic frequency lo = 10 (a), 
15 (6), 32 (c), and 38 rad/s (d) 



observed experimentally [||. In this case, the motion 
of the lower block is also periodic in time, but the 
time dependences of the parameters, which are given 
in Fig. 3, are not symmetric with respect to their 
zero values [3(|- The inset in Fig. 5,6 demonstrates 
the enlarged section conditionally marked by letter A, 
because this section has pronounced features, which 
cannot be distinguished in the main plot. Thus, the 
frequency u> affects the behavior of the tribological 
system in a non-trivial manner. By varying u), it is 
possible to select various modes of motion, which con- 
siderably differ from one another. Note that, at some 
frequencies, the stationary behavior of the system, 
which is established as a result of the system evolu- 
tion, depends on the initial conditions or the system 
prehistory. For instance, in Fig. 5,d, the initial value 
(po =/= gives rise to a mode similar to that exhib- 
ited in Fig. 5,c. This circumstance also confirms the 
presence of memory effects in the system, which were 
observed experimentally [jj. 

Figure 6 elucidates the influence of lubricant tem- 
perature X on the melting kinetics. The plotted de- 
pendences are divided into four sections. The tem- 
perature for the first section is the lowest, and, for 
every next section, the temperature increases, i.e. we 
have the inequalities Xi < X 2 < X 3 < X 4 . The de- 
pendence obtained in the first section, at X = Xi, 



reproduces the dependence shown in Fig. 3 in more 
details, because it was obtained at the same X-value. 
As the temperature is elevated to X = X2, the stick- 
slip mode of motion is realized, as it was at X = T\ . 
However, the maximum value of elastic stresses a e \ 
decreases at X = X2. As a result, the friction force F 
in the solid-like lubricant also decreases, as the tem- 
perature grows. As the temperature is elevated to 
X = X3, this tendency survives. Note that a reduc- 
tion of the sticking peak number with the tempera- 
ture growth is not a rule, and the opposite situation 
can take place. At X = X4, the lubricant is liquid- 
like all the time, and the elastic stresses equal zero. 
It is so, because, at this temperature, the condition 
X > X c ° (see Eq. (p|)) is obeyed even if e el = 0, i.e. 
the melted lubricant cannot solidify due to a reduc- 
tion of the relative shear velocity between the rub- 
bing surfaces. We do not know of any experiments 
devoted to similar researches of the influence of the 
temperature on the friction mode. Therefore, the de- 
pendences exhibited in Fig. 6 are a forecast. 

5. Numerical Experiment 

The dependences shown in Fig. 6 testify that the 
growth of the temperature X gives rise to a reduc- 
tion of the elastic stress amplitude a c \ and a reduc- 
tion of the friction force F maximum. Let is ana- 
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Fig. 6. Dependences of the coordinates X and x, the velocities 
V and v, the elastic stresses <r c i (see Eq. 113t ). and the friction 
force F (see Eq. I|10|l ) on the time t for the same parameters 
as in Fig. 3 and the temperatures Ti = 200 K, T 2 = 220 K, 
T3 = 250 K, and T4 = 300 K.The dashed curves correspond 
to the X(t) and V(t) dependences, and the solid curves to the 
x(t) and v(t) ones 

lyze the dependences of the <r e i and F amplitudes on 
the temperature T at various modes of functioning 
of the system in more details. We define the stress 
amplitude as cr c i,max := (cel.max - cr e i, min ) /2 and the 



friction force amplitude as F a 



/2, 



where 07 



and F m ,v are the maximum values of 



elastic stresses and friction force, respectively, and 
Cci.min and F mln arc their minimum values, which are 
determined within the complete period of parameter 
changes, T = 27r/w, after the stationary friction mode 
has been established. 

The dependences of the indicated quantities on the 
temperature are depicted in Fig. 7 for three types 
of lubricants: pseudoplastic (7 < 0), Newtonian 
(7 = 0), and dilatant (7 > 0) fluids. The upper panel 
demonstrates that, as the temperature increases, the 
elastic stresses er i,max decrease for all three types of 
fluids, i.e. the temperature elevation favors the lu- 
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Fig. 7. Dependences of the elastic stress, cr e i !max , and friction 
force, -Fmax, amplitudes on the temperature T for pseudoplastic 
(7 = —2/3), Newtonian (7 = 0), and dilatant (7 = 2/3) fluids 
as a lubricant. The parameters are the same as in Fig. 3 



bricant melting. Note that, for pseudoplastic fluids 
(7 = —2/3), which are used most often as lubricants 
in such systems, the stress amplitude attains maxi- 
mum values within almost the whole presented range 
of temperatures, but the melting occurs at lower T in 
this case. The lower panel of the figure shows the de- 
pendences of the friction force amplitudes -F max on the 
lubricant temperature T. It follows from the figure 
that the friction force decreases with the temperature 
growth only for pseudoplastic fluids, and it is mini- 
mal within the whole range of temperatures in com- 
parison with other types of fluids. For dilatant and 
Newtonian fluids and for the selected parameter val- 
ues, the maximum friction force does not change with 
the temperature growth. Since the elastic stresses for 
those fluids decrease as the temperature grows (the 
upper panel of the figure) , this means that the growth 
of T gives rise to an increase of the viscous compo- 
nent of the friction force, to which the second term in 
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fig. 9. Dependences of the elastic stress, cr c i max , and friction 
force, -F m ax, amplitudes on the spring stiffness constant K (see 
Eq. (0) for pseudoplastic (7 = —2/3), Newtonian (7 = 0), and 
dilatant (7 = 2/3) fluids as a lubricant. The parameters are 
the same as in Fig. 3 



formula (|10|) corresponds. In the situation concerned, 
this can happen only if the relative velocity of motion, 
V — v, increases. Note that, according to the figure, 
the amplitude of the friction force for the Newtonian 
and dilatant fluids remains constant when the tem- 
perature grows, even in the case a c \ = 0, i.e. when the 
friction force has only the viscous component. Since 
the F-amplitudes in the cases 7 = and 7 = 2/3 
coincide at all temperatures, it is not sufficient to ex- 
perimentally measure the total friction force in order 
to determine the friction mode. That is why the be- 
havior of the elastic, a e \, and viscous, cr v , stresses are 
additionally studied as a rule Note also that, 

according to the results demonstrated in Fig. 7, the 
application of pseudoplastic fluids is optimal to re- 
duce friction, because they favor the establishment of 
a mode with minimum force F, despite that the elas- 
tic stresses for such lubricants are maximum within 
almost the whole range of temperatures. 

To determine the dependence of the viscosity on the 
velocity gradient and the temperature, both real [28| 
and computer-assisted [29] experiments are carried 



out. The problem urgency is connected with the 
fact that the dependences of the viscosity on the in- 
dicated quantities are anomalous in the boundary 
friction mode in the case of nano-sized tribological 
systems. There can even be a mode, when the fric- 
tion force almost vanishes at cryogenic temperatures, 
which corresponds to a low viscosity of a lubricant 
and, accordingly, a very weak energy dissipation. In 
the English-language scientific literature, this mode 
was coined as "superlubricity" [4(|[4l|. Let us ex- 
amine the dependences of the friction force and the 
stresses for three types of lubricants - however, not 
on the temperature (as in Fig. 7), but on the propor- 
tionality coefficient k between the viscosity and the 
velocity gradient (see Eq. {7}). The corresponding 
plots are depicted in Fig. 8. Note that, in contrast 
to Fig. 7, different /c-values correspond to different 
lubricants, friction surfaces, or experimental geome- 
tries. This means that every point in the dependences 
exhibited in Fig. 8 corresponds to tribological sys- 
tems different by their properties. As one can see, for 
pseudoplastic fluids (7 = —2/3), the elastic stresses 
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remain constant with increase of the coefficient k. For 
Newtonian and dilatant fluids, the maximum stresses 
monotonously decrease with the increase of k. The 
friction force amplitude F max grows with the coeffi- 
cient k in the case of pseudoplastic fluid (7 = —2/3). 
At the same time, for the indices 7 = and 7 = 2/3, 
the friction force behaves identically as in Fig. 7, i.e. 
it remains constant. However, within the whole pre- 
sented range of k- values, F maK is minimal just for the 
pseudoplastic fluid; therefore, the latter is optimal for 
creating the conditions to reduce the friction in this 
case as well. 

In Fig. 9, the behavior of the examined quantities 
is illustrated, as the spring stiffness constant K in- 
creases. For the dilatant and Newtonian fluids, the 
elastic stresses <7 e i, ma x monotonously and slowly grow. 
In the case of the pseudoplastic fluid (7 = —2/3), the 
stresses drastically increase firstly, and afterward re- 
main almost constant. The friction force in this case 
(7 = —2/3) also grows to a certain value and, then, 
does not almost change. For the indices 7 = and 
7 = 2/3, the amplitudes of friction force F max lin- 
early increase with the spring stiffness constant K . 
and their magnitudes are equal as in the previous two 
figures. Hence, in this case, the pseudoplastic fluid 
also provides the minimum friction force in the sys- 
tem. Thus, a general conclusion can be drawn that 
the pseudoplastic fluids provide an optimal friction 
mode in the tribological system exhibited in Fig. 1, 
because the maximum friction force F max is the low- 
est for them. 

6. Conclusions 

In this work, a thermodynamic model was developed 
to describe the behavior of a tribological system func- 
tioning in the boundary friction mode. The model al- 
lowed a number of effects observed experimentally to 
be explained. It was shown that the stick-slip mode 
of motion is a result of the phase transition of the 
first order between the liquid- and solid-like states of 
a lubricant. The influence of the lubricant tempera- 
ture, the spring stiffness constant, and the coefficient 
of proportionality between the viscosity and the ve- 
locity gradient on the system behavior was analyzed. 
For pseudoplastic fluids, the elastic stresses and the 
friction force were found to decrease with the tem- 
perature growth. The increase of the spring stiffness 
constant induces the growth of the friction force and 
stresses for all types of lubricants. When the coef- 
ficient of proportionality k increases, the maximum 



stresses do not change substantially in the case of 
pseudoplastic fluids, whereas the friction force grows. 
For the sake of comparison, the results of calcula- 
tions obtained for the dilatant and Newtonian fluids 
were also reported. Modes, in which the displacement 
between the friction surfaces does not correspond to 
the direction of motion of the upper block, were re- 
vealed, which evidences the presence of memory ef- 
fects in the system. While developing the model, the 
thermodynamic potential with two stable stationary 
states was used, in which the zero and nonzero min- 
ima were separated by a maximum. However, it was 
found experimentally that the lubricant is character- 
ized by more than one type of transition and it can 
exist in a few (solid- or liquid-like) metastable states. 
For such a situation to be described, the additional 
terms of higher orders in the free energy expansion 
are sufficient to be taken into account. 
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ft. O. JlnmeHKO, A.M. 3acKona 

IJEPEPHBHACTHH PE>KHM ME>KOBOrO TEPT51 
5IK <5A30BHH IJEPEXm IIEPIIIOrO PO/IY 

P e 3 K3 m e 

Po3rji5myTO Tpii6ojiori t iHy CHCTeiviy, m,o CKJia^aeTBCH 3 ,z],box 

KOHTaKTyiOMIIX 6jIOKiB, O/THH 3 SKHX 3aKpilL7ieHHH M\yK A BO_ 

Ma npyjKHHaMH, a miniift npHBefleHnit b HenepepBHnit nepic-flii- 
hhhh pyx. /locjiifljKeHO KmeTHKy chctcmh b pe*ciiMi Me*coBO- 

TO TepTS, KOJIH M1>K aTOMapHO-rjiaflKHMH nOBepXHHMH 6jIOKiB 

3HaxoflHTBCH yjibTpaTOHKa njriBKa MacTHjia. /Lih onncy CTaHy 
MacTiijia 3anncaHO Biipa3 rjik rycTHHH bijibhoi eHeprii" y bh- 
rjijifli po3KjiaflaHH5i b psiA 3a CTeneHaMH napaivieTpa nopn^Ky, 
am 3BOflHTBCH flo MOflyjiH 3cyBy. IloKa3aHO, mp b miipoKO- 
My fliana30Hi napaineTpiB peani3yeTBCH nepepHBHacTnit pe>KHM 
pyxy, flo HKoro npiiBOflHTB nepioflHHHi <pa30Bi nepexc-flii nep- 
IHOrO pOfly Mi>K KiHeTHHHHMH pejKHMaMH Tepra. CnporH030- 
BaHO noBefliHKy CHCTeMH npn 3Mim 30BHiniHix Ta BHyTpinmix 
napaineTpiB. 
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